Electronic, dynamical and superconducting properties of CaBeSi 
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We report first-principles calculations on the normal and superconducting state of CaBeccSi2-cc 
(x = 1), in the framework of density functional theory for superconductors (SCDFT). CaBeSi is 
isostructural and isoelectronic to MgB2 and this makes possible a direct comparison of the electronic 
and vibrational properties and the electron-phonon interaction of the two materials. Despite the 
many similarities with MgB2 {e.g. a bands at the Fermi level and a larger Fermi surface nesting), 
according to our calculations CaBeSi has a very low critical temperature (Tc ^ 0.4 K, consistent 
with the experiment). CaBeSi exhibits a complex gap structure, with three gaps at Fermi level: 
besides the two a and tt gaps, present also in MgB2, the appearance of a third gap is related to 
the anisotropy of the Coulomb repulsion, acting in different way on the bonding and antibonding 
electronic tt states. 

PACS numbers: 74.25. Jb,74.25.Kc,74.70.Ad 



I. INTRODUCTION 

The complexity and the fragility of the effective pair- 
ing interaction - result of a fine interplay of opposite con- 
tributions - and its non-trivial dependence on chemical 
properties make the search of new superconducting mate- 
rials a very difficult task, even within the class of phonon 
mediated superconductors. As a result, research often 
goes along the lines of searching in the "neighborhood" 
of known superconductors. This has been the case of 
MgB2 that owes its superconducting (SC) propertie^i^i^ 
essentially to the presence of holes in the sp'^ a B-B 
covalent bonds, strongly coupled with the E2g stretch- 
ing mode. During the last years, MgB2 boosted experi- 
mental and theoretical research in the class of diborides, 
layered, graphite-like materials, all sharing some hope- 
fully relevant features with the "parent compound" . Al- 
though this effort did not succeed in finding supercon- 
ductors with better (or at least equivalent) properties 
than MgB2, it brought to the experimental discovery 
of some new superconductors like graphite intercalated 
compounds, boron doped diamond^ and to some the- 
oretical proposals 

Among MgB2-like materials, one interesting case is 
represented by CaSi2, which at ambient pressure has a 
rhombohedral structure and changes to a trigonal phase 
between ^ 7 and 9.5 GPa. At P > 16 GPa, CaSi2 has 
the AIB2 structure and SC critical temperatures (Tc) up 
to 14 K. The stability of the AIB2 phase at ambient pres- 
sure is achieved through hole doping. The stabilization 
process has been related^ to the filling of the tt antibond- 
ing bands (TTa). Let us consider a material CaXa;Si2-a; 
(X having a lower valence than Si): increasing x from 
zero introduces holes in the iVa bands, with a change in 
the X-Si bond that from an sp^-\ike character acquires a 
progressive sp'^ character, stabilizing the AIB2 structure 



Thus, hole-doped CaSi2 could represent (with some dif- 
ferences) a good candidate for an MgB2-like material, 
although its stabilization at high pressure poses some 
problems. Nevertheless, experimentally, the AIB2 phase 
was observed at ambient pressure by doping CaSi2 with 
^|io,ii Qj.^ ^ much less studied case, with BeJ^^ Re- 
cently, CaBexSi2-x was investigated in the doping range 
0.5 < X < 1^^ and was observed to have the AIB2 struc- 
ture at X = 0.75i2iii^ and, against the expectations, not 
to be a superconductor down to 4.2 

The goal of this paper is to clarify why this happens 
despite the presence of a holes, as in MgB2, and an even 
more pronounced nesting. We investigate the normal 
and SC phase of CaBeajSi2-a^ at x = 1 (CaBeSi in the 
following), in the ^4/^2 phase, with Be and Si atoms al- 
ternating within the honeycomb layers with an AA stack- 
ing. 

The paper is organized as follows: In Sees. [Ill and 
mil we sketch our computational framework and describe 
the computational details of our calculations. In Sec. 
HVlwe illustrate the normal state properties and the e-ph 
coupling of CaBeSi . The superconducting properties are 
discussed in Sec. [Vl Finally, in Sec. IVIi we summarize 
our results. 



II. DENSITY FUNCTIONAL THEORY FOR 
SUPERCONDUCTORS (SCDFT) 

The superconducting properties are studied within 
the density functional theory for superconductivity 
(SCDFT) ^ii^ a completely parameter free approach 
which allows to predict the SC gap and Tc values of real 
materialsJ^ii^ As the SCDFT method was discussed in 
full detail in previous papers ^li^ here we only sketch the 
main points concerning the theory, whose central result 
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is the generalized, self-consistent, gap equation: 




(1) 

where n and k are respectively the electronic band index 
and the Bloch wave vector; Ank is the SC gap function; (3 



is the inverse temperature; Enk = y {^nk — m)^ + \^nk\ 
are the SC quasiparticle energies, defined in terms of the 
gap function, the Kohn-Sham eigenenergies Enk of the 
normal state and the chemical potential /i. The universal 
kernel ICnk.n'k' appearing in Eq. ([T]) consists of two con- 
tributions /C = /C®~P^ + /C®~®, representing the effects of 
the e-ph and the e-e interactions respectively. /C^-p^ is 
temperature dependent and involves the e-ph coupling 
matrix elements \g'^\f and the phonon spectrum uOqjj^ 
while /C®~® contains the matrix elements of the screened 
Coulomb repulsion. Finally, the universal diagonal (and 
temperature dependent) term Znk plays a similar role 
as the renormalization term in the Eliashberg equations. 
We emphasize that Eq. ([1]) is not a mean-field equation 
(as in BCS theory), since it contains correlation effects 
via the SC exchange-correlation (xc) functional entering 
/C and Z. Furthermore, it has the form of a static equa- 
tion - i.e., it does not depend explicitly on the frequency - 
and therefore has a simpler structure (and computation- 
ally more manageable) than the Eliashberg equations. 
However, this certainly does not imply that retardation 
effects are absent from the theory. Once again, retarda- 
tion effects enter through the xc functional, as explained 
in Refs. dUi . 

An important feature of the SCDFT approach is the 
capability to include the n/c-resolved Coulomb repulsion 
V®~® ab-initio^ without any adjustable parameter. The 
Coulomb interaction was screened with a static dielec- 
tric matrix, within the Random Phase Approximation 
(RPA)ii(see below). 

III. COMPUTATIONAL DETAILS 

The electronic band structure Snk^ the phonon spec- 
trum cjqzy, and the electron-phonon (e-ph) and Coulomb 
matrix elements (ME) with respect to the Bloch 
functions necessary to solve the SC gap (Eq. [1]), 
were obtained within the planewave-pseudopotential 
methodic We used Troullier-Martins^ norm-conserving 
pseudopotentials,^- with 2s, 3s-3p and 3s-3p-4s for Be, Si 
and Ca, respectively, in the GGA-PBE^^ approximation 
for the xc functional. The electronic self-consistent cy- 
cle was performed with a 60 Ry energy cut-off, a 18^ 
Monkhorst-Pack /c-point mesh. 

The optimized lattice parameters (at x = 1 Be dop- 
ing) are a = 3.895 A and c/a — 1.112 , whereas the 
experimental values (at x = 0.75) are a = 3.94 A and 
c/a = 1.112^^ Our obtained, slightly smaller GGA- 
PBE a constant is justified by the larger amount of Be 
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c 


Si-Si 


Si-Be 


B-B 


CaBeSi 


3.895 


4.331 


3.895 


2.249 




MgBs" 


3.083 


3.52 






1.780 



^Experimental constants 
TABLE L Structural parameters (in A) of CaBeSi and MgB2. 



(with smaller covalent radius than Si) present in the cal- 
culated system. 

Phonons and e-ph ME were calculated via density func- 
tional perturbation theory. The phonons were com- 
puted on the irreducible set of a regular mesh of 8^ q- 
points, and a 16^ Monkhorst-Pack /c-points for electronic 
integration, with a smearing parameter of 0.35 eV. These 
parameters were sufficient to achieve convergence within 
0.5 meV in the frequency of the E2g{V) mode. Our cal- 
culated uoE^giX) = 59.0 meV compares well with the 
value u;E2g{^) = 57.1 meV obtained by frozen phonon 
calculations in Ref. 9, computed at the slightly larger 
a — 3.914 A. The e-ph ME were calculated on the same 
g-points grid as the phonons, and on a denser grid of 
24^ /c-points, while the RPA screened Coulomb ME were 
calculated^ on a mesh of 9^ x 9^ /c and k'-points. 
The SC gap function is extremely peaked around the 
Fermi surface (within the characteristic phonon energy), 
whereas at higher energies it is rather smooth (and it 
changes sign, due to the e-e interaction). This implies 
that a converged solution of Eq. ([T]) needs a denser Ap- 
points sampling around Ep^ while a coarser one is suf- 
ficient elsewhere. This highly non uniform mesh of the 
Brillouin zone is realized with 8 x 10^ and 500 indepen- 
dent /c-points for bands crossing and not crossing the 
Fermi level respectively. Finally, 15-20 self consistent it- 
erations were sufficient to achieve a complete convergence 
of the gap. 

IV. NORMAL STATE PROPERTIES 
A. Electronic structure and Fermi surface 

Fig. [1] shows the CaBeSi and MgB2 band structure. 
A general similarity is found: in both materials the <j, 
TT-bonding (tt^) and 7r-antibonding (tTq) bands cross the 
Fermi level {Ep)- As pointed out in Ref. 0, and sim- 
ilarly to CaAlSi^^ the reduced symmetry (space group 
P6m2) with respect to MgB2 (P6/mmm), related to the 
partially ionic nature of the B-Si bond, implies the split- 
tings of the cr and the Tr^-TTa bands at the K and H points 
of the Brillouin zone (BZ). In CaBeSi the tt^ bands are 
almost fully occupied, leaving only small holes pockets 
at K which give rise to the little tt spheres of the Fermi 
surface (FS) (see Fig. [2]). These replace the tt tubular 
structure present in MgB2. The TTa bands are only par- 
tially occupied, allowing the stabilization of the Si-Be sp^ 
network against a sp^-like distortion. In fact, the latter 
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FIG. 1: (color online) Band structure and DOS of CaBeSi 
(full lines) and MgB2 (dashed lines). Colors refer to different 
characters: a (blue), tt bonding (green), tt antibonding (red). 



FIG. 2: (color online) CaBeSi Fermi surface. The colors iden- 
tify the three gaps at Ef'- the largest a (yellow and blue 
cylindrical- like sheets); the intermediate tt bonding (green 
spheres); the lowest tt antibonding (violet). 



takes place in the presence of a larger amount of Si, i.e. 
with a larger filling of the TTa bands. ^ 

Due to the different electronegativity of Si and Be, we 
expect a change in the charge distribution related to their 
bond, in comparison with the B-B bond in MgB2. In 
Fig. [3] (upper panel) we plot the tt^ and TTa charge at the 
BZ H point. We see that the 7r5 (tt^) charge is clearly 
associated to the Si (Be) atom. Having in mind the dif- 
ferent occupation of the tt^ and TTa bands, we conclude 
that there is a charge disproportion in favor of Si. 

Another important chemical difference between MgB2 
and CaBeSi is the presence of Ca-d states at ~ 5 eV above 
Ep that strongly interact with tt and a bands, in different 
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FIG. 3: CaBeSi and MgB2 charge. Upper panel: CaBeSi 
TT bonding (left) and tt antibonding (right) charge at the H 
point; empty (full) circle refers to Si (Be). Middle and lower 
panels: CaBeSi (left) and MgB2 (right) a charge at the F 
point; in the middle (lower) panel, the z direction is perpen- 
dicular (parallel) to the page. 



ways along the zone. This contributes to the reduced 
in plane tt bandwidth (^ 1.1 eV and ~ 2.8 eV along 
FMK and ALH) in CaBeSi with respect to MgB2 5.5 
eV). In fact, while at the F point there is no interaction 
between tt states and the cation d orbitals, the interaction 
is possible in the high symmetry points A, M and 
therefore suppressing the large dispersion of these states 
observed in MgB2. A further reason for the reduced tt 
bandwidth in CaBeSi is a lower chemical pressure effect 
due to its larger unit cell in comparison to MgB2. 

Looking at the a charge in the two materials (Fig. [3l 
middle and lower panels), we notice a stronger localiza- 
tion in MgB2 along the B-B direction, whereas in CaBeSi 
the a charge has a clear ionic component and is more de- 
localized both in the in-plane and out-of-plane directions. 
This makes the Si-Be bond much weaker than the B-B 
one affecting the strength of the e-ph coupling. 

As far as the a bandwidth is concerned, the larger Si- 
Be distance (see Table H]) and a lower Si-Be interaction 
(compared to B-B) explain the reduced in plane a -band 
dispersion (for the lower a ~ 5 eV and ~ 8 eV in CaBeSi 
and MgB2, respectively). On the other hand, the larger 
out of plane dispersion of a bands in CaBeSi ( ^ 1 eV) 
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FIG. 4: (color online) Left panel: phonon structure of CaBeSi 
(full lines) and MgB2 (dashed lines). Right panel: CaBeSi 
phonon DOS (thin black line), Eliashberg function (thick 
brown line) and integration curve 2 F{ijj') / uj' duj' (dot- 
dashed green line) , whose final value is the total e-ph coupling 
A. 



versus MgB2 (~ 0.5 eV) is only partially justified with 
the larger 2;-extension of the CaBeSi <j charge (compen- 
sated by the ^ 25% larger interlayer distance), and is 
mainly related to the cr - Ca-(ia^2_^2 interaction allowed 
at the r point but not at A (see Fig. [T|). Difference in 
band dispersion gives rise to much of the warping of the 
corresponding Fermi surface. 



B. Phonons and electron-phonon coupling 

The previous discussion shows that, despite the gen- 
eral similarities, CaBeSi and MgB2 have rather different 
chemical and electronic properties. As expected, they de- 
termine both the dynamical properties and the electron- 
phonon coupling. In fact, CaBeSi frequencies are lower 
(see Fig.H]) in comparison with MgB2, mainly due to the 
larger mass of Ca and Si versus Mg and B. 

The mode is fairly flat along the in plane BZ sym- 
metry lines and it shows only a very weak renormaliza- 
tion along MF and AH lines (four times smaller than in 
MgB2)^ i^^i^^ due to the very small electron-phonon 
matrix elements. In turn, this is related to the delo- 
calized and ionic nature of the a bonds in CaBeSi (see 
Fig. [3|). In fact, the connection between strongly cova- 
lent bonds and strong e-ph coupling seems to be a general 
featurei^^^ 

The Bi^ mode lower than the E2^ everywhere in the 
BZ, without exhibiting the features found in CaAlSi, 
where it is very soft, due to a strong interband coupling 
between the interlayer and TTa states 

The strongly reduced e-ph renormalization of the 
CaBeSi E2^ mode is not related to poor FS nesting fea- 
tures in this material, but only to the small value of the 
e-ph ME themselves. In fact, we have calculated the 
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FIG. 5: (color online) CaBeSi superconducting gap as a func- 
tion of the energy distance from the Fermi level. 



cr nesting function ry^ = j^^Y^k^a ^i^k) S{ek+q) (where 
k -\- q ^ a and Na{0) is the a DOS at Ef), obtaining larger 
values for CaBeSi than for MgB2 (roughly a factor of 2). 

The calculated total e-ph is A = 0.38, which makes CaBeSi 
a weak coupling superconductor, comparable to Al or Mo, 
with a Tc = 0.4 K. The two-band resolved values (see Eq. 
El) are X^a = 0.29; A^^ = 0.21; A^,, = 0.15; A^^ = 0.12 
whereas the MgB2 values are^° Xaa = 0.83; Xair = 0.22; X^ra — 
0.16; Atttt = 0.28. The comparison further shows the dramatic 
reduction of the g-g coupling in CaBeSi. 

Although weaker than in MgB2, the contribution from the 
E2^ mode is important also in the CaBeSi Eliashberg function 
(Fig. 2]), where we see a peak at ~ 60 meV which strongly 
enhances a corresponding structure in the phonon density of 
states. This peak gives roughly one fourth of the total e-ph 
coupUng, while the main contribution to the coupling comes 
from low frequency modes. 

We note that charge localization not only enhances the e- 
ph interaction, but also the Coulomb ME,— which are lower 
in CaBeSi than in MgB2. However, this reduction is not so 
dramatic as for the e-ph ME. This is due to the different 
structure of the ME of the two interactions (see discussion 
below) . 



V. SUPERCONDUCTING STATE 

The solution of the self-consistent gap equation (Eq. ([T])) 
including the anisotropic e-ph ME 1^'^^/ and, the nk- 
resolved RPA Coulomb matrix elements VV^/w, reveals an 
unexpected complex structure with clearly separated three 
gaps at Ef (Fig. [5)). The calculated critical temperature is 
very low (Tc= 0.4), lower than the upper limit (4.2 K) set by 
the experimental results. Unlike in MgB2, in which super- 
conductivity is interpreted within a two-band model,— i^*^*^ 
in CaBeSi there is a further TVb-TTa gap splitting. As in MgB2 , 
the largest gap is related to the a FS sheets (cylindrical-like 
structures in Fig. [2}, the intermediate one to TVb sheets (small 
hole spheres) and the lowest to TVa sheets. The additional ivb- 
TVa gap splitting is a peculiar feature of CaBeSi not present in 
MgB2 in which the two a and the two tt gaps merge together. 
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In order to understand the origin of this sphtting we perform 
some additional computational experiments, solving the gap 
equation (i) completely neglecting the Coulomb interaction, 
(ii) including only the averaged Coulomb term 

N(^e)N{d) ^ ^nk^n'k^ X ^(^nk - c) ^(e^/fc/ - e) . (2) 



and (Hi) with isotropically averaged Coulomb and phononic 
interactions, corresponding to the dirty limit. 

In both (i) and (ii) cases, the three gap structure is de- 
stroyed, bringing back to a two-band, MgB2-like gap struc- 
ture. In case (in), instead, superconductivity is completely 
lost. As a result, we predict superconductivity in CaBeSi only 
if the anisotropic structure of the interactions is included. In 
the real system, very likely, disorder on the Si-Be sublattice 
can produce interband Tv-a impurity scattering, therefore re- 
ducing Tc further. 

In the following we will analyse the three gap structure, 
with particular emphasis in understanding the TTb-TTa splitting 
observed. To this purpose, we perform a four-band model 
analysis, splitting the Fermi surface in: internal a band (cri), 
external a band (cr2), TVb and TVa bands. 

We have calculated (see Table Ulla and b ) the corresponding 
density of states Nn and the BCS-like e-ph couplings, X^n' 
and y""^^, where: 



a dos 




cr2 


7Tb 


7Ta 


N 


0.12 


0.38 


0.1 


0.58 




b e-ph 






7Tb 


7Ta 


cri 


0.07 

y\j.OZ ) 

0.07 

(0.61) 


0.23 

0.22 

(0.57) 


0.03 

yyj.Zo ) 

0.02 

(0.23) 


0.20 

0.19 

(0.32) 


7Tb 


0.03 

(0.28) 


0.09 

(0.23) 


0.03 

(0.35) 


0.16 

(0.28) 


TTa 


0.04 

(0.34) 


0.12 

(0.32) 


0.03 

(0.28) 


0.07 

(0.12) 



V 



e — ph 



(3) 



c el-el 


cri 


(72 


7Tb 


7Ta 


cr2 


0.11 

(0.91) 
0.08 

(0.71) 


0.27 

(0.71) 
0.32 

(0.84) 


0.03 

(0.30) 
0.03 

(0.29) 


0.15 

(0.26) 
0.13 

(0.23) 


7Tb 


0.04 

(0.30) 


0.11 

(0.29) 


0.09 

(0.94) 


0.20 

(0.35) 


7Ta 


0.03 

(0.26) 


0.09 

(0.23) 


0.03 

(0.35) 


0.51 

(0.86) 



and TABLE II: Table a: Band resolved density of states 

[st/(cell-eV)]. Table b: Band resolved e-ph couplings A^^/ 
I nn' |2 (in bold) and BCS effective potentials V^^~/^^ = Xnn' /N'^ (in 

ye-Tph _ ^ \9k,k' parcuthcsis [eV]). Table c: Band resolved el-el interactions at 

NnNn' j^^^j^,^^, ^ ^k'-k,u the Fermi energy /Jnn' (in bold), and el-el matrix elements 



X S{£r^k - EF)5{er^'k' " Ep). (4) 



^nn'^ at the Fermi energy (in parenthesis [eV]). 



First, we notice that the a submatrix V^~^ is very homoge- 
neous. Second, the g-tt scattering gives the same contribution 
to both cr gaps, which are then identical. Superconductivity 
in the tt states is more complicated, essentially because, un- 
like <7 bands, TTa and iTb bands originate from different sublat- 
tices. As a consequence, the tt submatrix is not homogeneous. 
Moreover, the tt-g interaction is of the same size of the tt-tt 
one, and TTa and TTb have different density of states, being 
A^TTa — 6 • A^TTb • order to understand the qualitative struc- 
ture of the SCDFT results, we considered a BCS-type approx- 
imation of the 4-band Eliashberg equations.^^ This model cal- 
culation confirms that the inclusion of the average Coulomb 
interaction (Eq. [2]), does not produces TTb-TTa gap splitting. 
In fact, the splitting is recovered only considering the (band) 
anisotropic V^~,^ Coulomb ME, reported in Table Hlb. These 
are, as for the el-ph coupling, roughly homogeneous in the a 
submatrix and therefore not able to split the a gap, but in 
the V^~^ submatrix, the intraband interaction is ~ 2.7 times 
larger than the interband interaction (that couples states in 
different sublattices) . The huge difference between V^~t^^ and 
V-^a-^a (that makes the /iTTbTr;, and 

ligible), lead to a much stronger suppression of A^r^ relative 
to Attj, , eventually bringing to the three gap structure. 



VI. SUMMARY AND CONCLUSIONS 

We have calculated the normal and superconducting state 
properties of CaBea;Si2-cc {x = 1) in the AIB2 phase, within 
the density functional theory for superconductors (SCDFT). 
The chosen doping level is not far from the experimental dop- 
ing (x = 0.75) where this phase is found stable and homoge- 
neous. CaBeSi is isostructural and isoelectronic to MgB2 and 
the electronic, vibrational properties and electron-phonon in- 
teraction of the two materials are compared directly: While 
the band structures present strong similarities, with both a 
and TT bands crossing the Fermi level, the phonon structure 
and the e-ph interaction differ substantially. In particular, the 
less localized a charge of CaBeSi brings to a dramatic reduc- 
tion in the E2g electron-phonon coupling, with a consequent 
reduction of the phonon renormalization seen in MgB2. This 
fact makes CaBeSi a weak coupling superconductor with e-ph 
A = 0.38 and Tc = 0.4 K, in spite of a nesting at the Fermi sur- 
face twice as big as in MgB2 . Interestingly, CaBeSi exhibits 
three superconducting gaps at the Fermi level. While, as in 
MgB2, the a-TT gap splitting is related to the different e-ph 
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coupling in these bands, the further nb-TVa sphtting is a pure 
effect of the complex structure of the anisotropic Coulomb 
repulsion, acting in different way on the TVb and TVa states. 
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